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ABSTRACT 

We have discovered a large number of circular and elliptical shells at 24//m around luminous central 
sources with the MIPS instrument on-board the Spitzer Space Telescope. Our archival follow-up effort 
has revealed 90% of these circumstellar shells to be previously unknown. The majority of the shells is 
only visible at 24/im, but many of the central stars are detected at multiple wavelengths from the mid- 
to the near-IR regime. The general lack of optical counterparts, however, indicates that these sources 
represent a population of highly obscured objects. We obtained optical and near-IR spectroscopic 
observations of the central stars and find most of these objects to be massive stars. In particular, we 
identify a large population of sources that we argue represents a narrow evolutionary phase, closely 
related or identical to the LBV stage of massive stellar evolution. 



Subject headings: ISM: bubbles — stars: circumstellar matter — stars: early-type 
line. Be — stars: mass loss — stars: Wolf-Rayet 



stars: emission- 



1. INTRODUCTION 

The Spitzer Space Telescope has provided us with an 
unprecedented view of the mid-IR sky in terms of both 
resolution and sensitivity. Invariably, the opening of ev- 
ery new wavelength window carries with it the potential 
for the discovery of previously unrecognized phenomena 
and source populations, which have gone unnoticed at 
other wavelengths. Motivated by our discovery of an un- 
usual well-defined elliptical shell - lacking a central en- 
ergizing source and only visible at 2 4^(m - in some of th e 
first data from the Spitzer mission ()Morris et al.ll2006l ). 
we executed a search for similar objects that might shed 
light on its still puzzling origin and nature. 

We utilized the publicly available 24/xm data products 

i delivery Version 3 .0) of the MIPSGAL Legacy project 
Carev et"ani2009D obtained with t he Multiband Imag - 
ing Photometer for Spitzer (MIPS. iRiecke et al.|[200l . 
We conducted our search by eye, inspecting all the im- 
ages of the inner 248 square degrees of the Galactic plane 
covered by the MIPSGAL survey. Although there was no 
prediction of whether 24/im shells would be rare or com- 
mon occurrences, we were still surprised by our discovery 
of a large number of prominent, more or less ring-shaped 
24/im objects: we find over 200 such shells in the MIPS- 
GAL survey area. Unlike the shell that originally mo- 
tivated the search, however, the majority of the objects 



wachter@ipac.caltech.edu, mauerhan@ipac.caltech.edu, vandykii 



skatka@dtm.ciw.cdu 
pmorris@ipac.caltech.edu| 

Visiting Astronomer, Cerro Tololo Inter- American Observa- 
tory, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. 



have some indication of a central source, at times quite 
luminous even at 24/(m. 

Follow-up investigation utilizing the complementary 
data from the GL IMPSE Spitzer Legacy project 
(|Beniamin et al.|[2003f) which covers the same survey area 
with the IRAC 3.5-8.0/im bands, as weh as 2MASS 
and the Digitized Sky Survey (DSS), shows that about 
60% of these shells are only detected at 24/(m (see Fig- 
ure [T]). This is somewhat unusual, as we generally expect 
a strong emission component due to polycyclic aromatic 
hydrocarbons (PAHs) at 8/im if we are simply observ- 
ing warm dust co ntinuum emis s ion. In the case of the 
shell discussed by iMorris et al.l (|2006D . the lack of 8/tm 
detection is explained by a pure line emission spectrum 
with a dominant [OIV] 25.89/tm line. We speculate that 
some of our 24/tm-only shells might exhibit similar mid- 
IR spectral characteristics. 

Ring-like shell morphology is observed frequently in 
planetary nebulae or in material associated with a vari- 
ety of massive stars, such as AVolf-R ayet (WR) stars and 
red and blue supergiants (jChulll99TI ). So far, most known 
circumstellar rings and shells have largely been studied 
using narrow band optical filters (e.g. Ha, [OIII]) and, 
despite being labeled as ring-shaped, many of these nebu- 
lae appear quite amorphous. For our Spitzer-discovered 
shell sources, a SIMBAD search within 2' of the shell 
locations revealed that most (90% of the shells, 80% of 

)ipac.c«ifecSfifeaSf;h8&^^ ^^'^^ ^'^^ ^^^^ previ- 

ously sLudied. Vve presenl here the results of our efforts 

to determine the nature of the central stars that produce 
the 24/im shells and to characterize the properties of this 
substantial population of new sources. 
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2. SOURCE SELECTION 

Identifying a region of 24/im emission as a more or less 
symmetrical shell is necessarily a subjective undertak- 
ing. In order to compile the most homogeneous sample 
possible we have selected a subsample of our shells for 
dedicated follow-up observations, based on the following 
criteria: 

• The 24/im shell has a symmetric (circular or ellip- 
tical) shape with a well defined boundary. 

• The presence of a likely central source, detected at 
a minimum of two different wavelengths among the 
24/im, 8.0/im, 3.6/im, 2MASS Ks and J images. 

The second criterion is particularly tricky, since we are 
dealing with the crowded fields in the Galactic plane and, 
hence, a relatively high a priori probability for chance 
superpositions. In order to classify a given source as a 
"likely" central source, we first fitted each shell by eye 
with a circle or ellipse, then searched for sources within 
1" of the geometric center of the fitted shape. Shells 
with bright 24/im central sources represent the most re- 
liable cases, since the overall source density at this wave- 
length is significantly lower than in the shorter wave- 
length bands. In other cases the central source identi- 
fication is more tentative, and we comment on this for 
individual sources in Section 8. Larger shells have greater 
uncertainties in the determination of the geometric cen- 
ter purely from the error inherent in the fitting process. 
In addition, the literature certainly contains examples 
of shells excited by sources that are significantly offset 
from the geometrical cent er (e.g. [Whitehead et al.lll988l : 
iSiebenmorgen et al.|[l998l ). 

Our list of targets is presented in Table [H For each 
shell we have indicated the 2MASS designation of the 
likely central source, the approximate angular size of the 
shell, the presence of the central source counterpart at 
various wavelengths, and whether the 24/xm shell is also 
detected at 8/im. We searched for known objects at the 
position of the central sources in SIMBAD, using a 2-5' 
radius depending on the size of each shell. The results 
are listed in the last column of Table [U as well as the 
results of our IR and optical spectroscopic follow-up of 
the central sources. 

In order to facilitate the correlation of central source 
properties with the morphological characteristics of the 
shells, Spitzer 24 and 8/im, 2MASS Ks and DSS R band 
images of the shells are displayed in Figure [T] (in its en- 
tirety available only in electronic format). The assumed 
central source is marked. 

3. OBSERVATIONS OF THE CENTRAL SOURCES 

While some of the central sources from Table [T] can 
be matched to objects in SIMBAD, the vast majority 
has not been observed before. It is impossible to deter- 
mine the nature of the central sources based solely on 
archival photometric data, since the extinction ~ both 
intrinsic to the source and between us and the source - 
and the distance to the shells are unknown. Hence we ob- 
tained optical and IR spectroscopy of the central sources 
of these newly discovered shells to identify their nature. 
The spectroscopic data obtained is indicated by "IR" for 
infrared and "OPT" for optical in Table [l] respectively. 



3.1. Near- Infrared Spectroscopy 

Near-Infrared spectroscopy of a number of the central 
stars of our shells was obtained on 2009 June 23-25 UT 
with the Ohio State Infrare d Imager and Spectrograph 
(OSIRIS; iDepov et al1ll993D on the Southern Observa- 
tory for Astrophysical Research (SOAR) 4.1m telescope. 
Our low- resolution spectra (R~1200) were obtained us- 
ing the 1" cross-dispersed slit and the fS camera which 
covers the J, H and K bands simultaneously. Spectra 
were acquired in an "ABBA" nodding sequence for sky 
subtraction. We utilized sky spectra for wavelength cal- 
ibration. Data reduction and spectral extraction and 
calibration was performed in the standard manner us- 
ing IRAF. For telluric correction, spectra of bright AOV 
standard stars ( Ks ~ 6-7 mag) wer e obtained. The IDL 
program xtellcor (jVacca et al.ll2003l ) was used, which re- 
moves model H i absorption lines from the AOV spectra 
before application to the science data. Observing condi- 
tions were photometric on the first night, but strongly af- 
fected by variable clouds during the following two nights, 
which limited the observations to the brighter targets in 
our source list. Note that no J band spectra could be 
extracted for stars #4, 5, 6a, 6b, 8a, 14, 20, 22, 29, and 
44. 

The IR spectrum for star #52 was obtained with 
the Palomar Hale 200" telescope and the TripleSpec 
(TSPEC) near-infrared spectrograph on 2009 July 14 
UT. TSPEC provides simultaneous coverage of the J, H, 
and K bands (A = 1.0-2.4 /xm) and produces a moderate 
resolution spectrum [R k. 2500-2700) through a 1" x 30" 
slit. Spectra were obtained utilizing an ABBA telescope 
nodding sequence for sky subtraction and bad-pixel sup- 
pression. Wavelength calibration was performed using 
the OH emission lines in the sky spectra. The AOV star 
HR 6958 was observed as a standard, and a telluric spec- 
trum was applied to the science data using the program 
xtellcor. 

3.2. Optical Spectroscopy 

For the limited number of sources with optical counter- 
parts, optical spectroscopy was obtained with the 200" 
Hale telescope at Palomar Observatory on 2008 Septem- 
ber 3 and 4 UT using the Double Spectrograph. This 
low to medium resolution grating spectrograph uses a 
dichroic to split the light into separate red and blue chan- 
nels which are observed simultaneously. The blue side 
CCD has 15/im pixels and a scale of 0.389"pixel~^, the 
red side CCD 24/im pixels with a scale of 0.468"pixel~'^. 
We used the 316 lines/mm grating in first order, dichroic 
D55 and a slit width of 1". The resulting spectra cover 
4000-5600A and 5800-8300A with a dispersion of 2.0 and 
2.4Apixel~^ on the blue and red side respectively, and 
have a resolution of 5-7A. For wavelength calibration, 
FeAr (on the blue side) and HeNeAr (on the red side) 
arc lamp spectra were obtained at each new telescope 
position. Bias and flat-field corrections to the raw spec- 
trum images were performed with IRAF. Various rou- 
tines within the IRAF packages twodspec and onedspec 
were used to extract and wavelength/flux calibrate the 
spectra. Both nights were photometric with variable see- 
ing. Observations of spectrophotometric standard stars 
were carried out, however, these were mainly used to 
remove the instrumental response function. A signal-to- 
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noise ratio of 40-100 was achieved in the final spectra. 

4. SPECTRAL CLASSIFICATION 

4.1. IR Spectra 

At first glance, our central source IR spectra can be 
roughly divided into two main groups: those that show 
emission lines or early spectral type features and those 
that exhibit absorption features indicative of late type 
stars. 

4.1.1. Late Type Stars 

For the late type stars, we determined the spectral 
type and luminosity class based on the strength of the 
CO bandhea,d at 2.29/z m, following t he me thodology of 
iFiger et all ()2006[ ) and iDavies et al.l (|2007| ). Both ways 
provided consistent results. We also compared the ap- 
peara nce of the spectra t o those from the IRTF spectral 
atlas (jRavner et al.ll2009l ). In addition, we performed the 
same analysis for selected spectra from the atlas (of stars 
with known spectral type and luminosity class) as an in- 
dependent check for the validity of our measurements. 
We estimate that our derived spectral types are accurate 
to within ±1 subtype. The results are listed in Table [2] 
and a selection of our late type spectra are displayed in 
Figured] Note that several sources show Br7 in emission. 

We caution that veiling of the absorption features by 
a hot dust continuum could affect our spectral type de- 
termination. This would lead to an underestimate of our 
measured equivalent width, and hence the true spectral 
types of our sources would be later than those we derived. 
In order to investigate the presence of circumstellar dust, 
we plotted the 2MASS colors of the stars from Table [2] 
in a near-IR color-color diagram (not shown) . Typically, 
IR excess due to circumstellar material will cause a shift 
to redder H — Ks colors relative to the unreddened stel- 
lar loci. We find that all of our stars fall on the red 
(super)giant reddening line within the photometric un- 
certainties. The most deviant point corresponds to star 
#23b, which is unlikely to produce the associated shell 
as discussed below. 

As indicated in Table [21 each CO bandhead equivalent 
width measurement corresponds to a solution for both 
giants and supergiants, and we are unable to break the 
degeneracy solely based on our spectra. Hence we in- 
vestigated whether photometric information could pro- 
vide additional constraints to distinguish between these 
possibilities. For each source, we first calculated the 
total (interstellar and circumstellar) Ay by comparing 
its 2MASS photometry (provided in Table [3]) with the 
unreddene d J — K colors of the appropriate spectral 
type from iCoxl (1200 0). utilizing the reddening law by 
ICardelh. Clavton. fc Mathisl (|1989( ). Finally, we derived 
the distances and the physical sizes of the surrounding 
shells for the two possible luminosity classes, respectively. 
These results are also listed in Tabled 

For all stars, the derived Ay exceeds the average 1 mag 
kpc~^ value by several factors, implying that the line of 
sight to these sources intersects localized areas of high 
extinction. Minimizing this excess extinction would in- 
dicate the supergiant case for all sources. However, for 
star #26b we favor the giant scenario, since the distance 
of >23 kpc for the supergiant case results in an average 
extinction of only 0.4 mag kpc~^ and the largest physical 



shell size of the sample. We derive and list these param- 
eters for all of the sources for completeness, however, 
stars #8b, #23b, and #28b are probably not producing 
their respective observed shells. These particular stars 
are members of apparent doubles, where the center of the 
fitted shell geometry coincides with the other component 
of the double. While we cannot exclude the possibility 
that the off-center source is responsible for the shell, it 
appears more likely that the central source is producing 
the shell, in particular in the case of #23b and #28b, 
where the central source was found to be a hot star (see 
the following section). 

Most of our knowledge about large scale circumstel- 
lar shells around red giants a nd supergiants dates back 
to the IRAS 6Q/xm s urveys of [St encel et aj] (|1989[ ) and 
lYoung et all ()1993l ). lYoung et a l. (1993) surveyed 512 
red giant stars, 15% of which were found to have shells 
with radii ranging from 0.1-4.6 pc. The average shell 
radius was 0.74 pc, with radii of > 2 pc being quite rare. 
iStencel et al.l ()1989D conducted a similar survey for 111 
red supergiants, 25% of which showed evidence for shells 
with radii between 0.1-4 pc (most clustered between 0.1- 
1.5 pc). Our shells have very similar size characteristics, 
which is somewhat surprising given that, a) our shells 
were observed at 24/xm and one might expect size differ- 
ences between the shells at various wavelengths, and, b) 
the resolution of the IRAS 60/^m images was quite limit- 
ing, such that these surveys were biased towards detect- 
ing only the largest shells of any underlying distribution. 

4.1.2. Emission Lme / Early Type Stars 

Most of our IR spectra with emission lines resemble 
those of the so-called "transition objects" , massive stars 
of type Of, WNL, Be, B[e] and Luminous Blue Variable 
(LBV), thought to span the evo lutionary stages be tween 
O and Wolf-Rayet (WR) stars (jMorris et al .lll99g). The 
classification of these stars based on IR spectra is noto- 
riously difficult. We find that our spectra separate fairly 
cleanly into two distinct groups simply based on the ap- 
pearance of their spectral features. When considering 
commonalities, we give less weight to the appearance of 
the He i 2.058/im line, since it is highly variable due to 
its susceptibility to optical depth effects. 

The first group is comprised of stars #1, 3, 10, 13, 
36, and 52. For star #36 we only have a K band spec- 
trum, J and H were too noisy to reliably identify fea- 
tures. All stars exhibit strong He i 1.70 and 2.112/im, as 
well as prominent H i emission in all three bands (Fig- 
ure [3]). For all figures, t he line iden t ifications have bee n 
com piled from Morris efal] (|1996D . iFiger et al.] ()1997D . 
and lCrowther et al.i ()2006l ). We have divided the second 
group into two subgroups, distinguished by the strength 
of the H I features. The spectra in group 2A (stars #23a, 
24, 45, and 46) are still dominated by H i emission of sim- 
ilar strength as group 1 in the J, H, and K band, but 
the striking difference with the spectra of group 1 lies in 
the strong Fe ii and Mg ii emission and weak or absent 
He I (1.70, 2.112^m) features (Figure H]). In addition, 
[Fe ii] is detected in all the stars except #45. Group 
2B (stars #lla, 14, 17, 29, 32, 44) is characterized by 
much weaker H i emission (or even weak H i absorption) 
compared to the first two groups. However, similar to 
the stars from group 2A, the spectra exhibit Fe ii or Mg 
II emission lines (Figure [5]). We regard group 2B as an 
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extension of group 2A towards weaker H i features. In 
particular, star #29 could have been assigned to either 
group 2 A or 2B, and emphasizes the connection between 
the two groups. Star #32 is included in group 2B based 
on the similarities of its spectral features in the H band 
and the presence of Pa/3 emission in the J band; how- 
ever, it has no obvious emission lines of Fe ii, and the K 
band is fe atureless. In SIMBAD , the source is classified 
as M2: bv lRaharto et all (|1984D . but an M spectral type 
can be ruled out based on our IR spectrum, given the H 

1 absorption spectrum in the H band and the lack of CO 
absorption features in the K band. 

We inspected all of our spectra for the presence of the 
He II 2.189/im line, since this is one of the key lines in 
identifying WR and O type stars. Only star #10, #52, 
and #llb (Figure [6|) exhibit convincing He ii 2.189/im 
emission. Here, care must be taken to not mistake He i 
2.185/im emission for that of He ii 2.189/im. Our spectra 
have sufficient resolution to clearly distinguish between 
the two lines, e.g., star #1 shows He i 2.185 and not 
He II 2.189/im emission. The ratios of the equivalent 
widths (EW) of certain emission lines in the i^-band are 
commonly used a s a WR subtype diag nostic (e.g., see 
iFiger et all 119971 iCrowther etall 120061 ). For star #52 
and #10 we measure EW(A2.189 /im)/EW(A2.165 /im)= 
0.5 and < 0.1, respectively, which correspond to WR sub- 
types of WN7 and WN9h. We also have optical spectra 
of star #52 which independently confirm its WN7 clas- 
sification based on a different set of spectral criteria (see 
Section [421). For star #llb we find EW (Civ 2.076) / 
EW (Cm 2.110) < 1, consistent with a WC9 subtype. 

Unfortunately, no classification scheme based on the 
EW measurement of features exists for the remaining 
spectra, and we have to re ly on compar i son to the 
various spectral at lases, e.g., iHanson et al.l ()1996[ ) and 
iMorris et al.l (jl996( ). to classify these sources. The stars 
from our group 2 strongl y resemble th e spec tra of LBV 
and Bc/B[e] stars from iMorris et al.l (|1996f) shown in 
their Figures 5, 6, and 7. The exception in that sam- 
ple is the LBV He 3-519, which does not exhibit Fe ii 
emission, but instead looks like the spectra of stars in 
our group 1. It is interesting to note that [Fe ii] 1.644/im 
emission has been found to be a common property of the 
shell emission among LBVs with prominent nebulae, but 
has rarely be en seen in the spectra of their central sources 
(jSmithll2002l ). In contrast, we observe [Fe ii] emission in 
about half the sources of group 2. Based on these char- 
acteristics, we tentatively classify the stars in our group 

2 as Be/B[e]/LBV candidates. This is supported by the 
fact that star #29, #44, and #45 are identified as LBVs 
or LBV candidates in the literature (for more details see 
Section [4X3]) . 

Among the stars of group 1, we have already classified 
source #52 as WN7 and #10 as a WN9h star. The 
K band spectra of the remaining four stars in group 
1 (#1, 3, 13, 36) closely match the appearance of the 
WN9 source, except for the presence of the He ii 2.189/im 
line. Several studies have classified similar-looking spec- 
tra as WN9, despite the lack of He ii 2.189/tm emission. 
For example, iFiger et all (|1999f ) show several spectra of 
Quintuplet cluster stars that are virtually "clones" of our 
group 1 spe ctra, and none exh ibit evidence for features 
at 2.189/im. iShara et al.l ()2009[ ) display a K band spec- 
trum of the known WR 102d, classified as WN9, that 



equally lacks He ii 2.189/im. Finally, our spectra at 
both H and K ba nd resemble the Ofpe/WN9 stars of 
IMorris et all (|i996') (th eir Figures 2 and 3). Only some 
of those s tars ex hibit He ii 2.1 89jLtm. However, the Of 
spectra of IHanson et all ()1996n demonstrate that He ii 
2.189/im is always present for spect ral types earlier than 
09. Incidentally, the Oe spectra of IHanson et al.l ()1996f ) 
are very similar in appearance to our group 1 stars and 
do not show He ii 2.189/im. While the exact classification 
of these types of spectra is clearly still an open question, 
they undoubtedly indicate massive stars most likely on 
their way to a WR phase. For the purposes of our study, 
we will label these stars as Oe/WN9. 

The final two spectra remaining to be classified are 
shown in Figure [6l together with our newly discovered 
WC9 star. Star #20 exhibits Br7 emission. He i 2.112/tm 
in absorption, as well as possible Mg ii emission in the K 
band and He i 1.70/im absorption and a hint of the H i 
series in absorption. Based on the absence of obvious He 
II features, we classify star #20 as type Be. Source #28a 
also shows He 1 1.70 and 2.112/im (and possibly 2.185/im) 
in absorption and shallow absorption of H i. Based on 
the K band spectrum and comparison to IHanson et al.l 
(1996), the spectrum of star #28a most closely resembles 
that of HD 168021 (BO lb), HD 213087 (BO. 5 lb) and HD 
226868 (HMXB, 09.7 lab) in the relative strength of the 
features at Br 7 and 2.112/im. However, in the H band 
the shallow and broad H i features indicate luminosity 
class V and, together with the str ength of He i 1 .70/im 
line, a spectral type of O9-B0.5 ([Hanson et al.l fl 9981. 
We classify star #28a as OB, noting that the luminosity 
class issue still remains to be resolved. The final spectral 
classification for each source is listed in Table [TJ 

4.1.3. Known LBVs 

Sources #29 (WRAY 17-96 = Hen 3-1453), #44 
(G024.73-h00.69 = V481 Set), and #45 (G026.47-h00.02) 
are identified in SIMBAD as a LBV candidat e 
(lEga n et al.l[200l) . confirmed LBV (IClark et al.l [20051). 
and LBV candidate (IClark et al.ll2003f ). respectively. All 
have been found to be surrounded b y ring nebulae in t he 
mid-IR with ISO (Egan et al. 2061IClark et^^M^ . 

LBVs are generally defined as a class of massive, unsta- 
ble stars that populate the uppe r left corner of the HR di- 
agram (for a detailed review see lHumphrevs k, DavidsonI 
119941 ). These high luminosity stars are thought to be 
a short-lived phase before the WR stage, where a mas- 
sive star sheds a significant amount of its mass. Conse- 
quently, LBVs are characterized by high mass loss rates, 
significant photometric and spectral variability and oc- 
casional mass eruptions, the most famous example be- 
ing the historic 10-magnitude brightening of rj Carinae 
in the 1840s. These outbursts are extremely energetic; 
they have been linked to the discovery of rapidly ex- 
panding nebulae surrounding LBVs and, in several cases, 
have pre viously been mistaken for supernovae in external 
galaxies ([Van Dvk et al.ll2000D . More common than such 
giant eruptions, however, is the presence of a circumstel- 
lar nebula or shell, indicating more modest, past mass 
ejection episodes. Out of the 12 confirmed LBVs, seven 
are associated with circumstell ar ejecta, as well as 13 of 
the 23 LBV candidates listed bv lClark et all ([2005D . The 
underlying physical cause that triggers the LBV instabil- 
ity is not yet understood, the current consensus invokes 
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instabilities in the outer layers of the stars as they ap- 
proach the Eddington luminosity during their evolution. 

Our newly acquired spectra of the known LBVs allow 
us to search for the expected spectral variations in com- 
parison to the spectra previously published in the litera- 
ture. LBVs in quiescence (at visual minimum) typically 
exhibit spectral types of Be or B[e], rich with emission 
lines of hydrogen, helium, iron, and sodium. In many 
LBVs, the Fe ii and [Fe ii] emission lines are most pro- 
nounced during the quiescent state. During an outburst, 
the spectra may morph into those of much cooler type 
A/F supergiants as a pseudo-photosphere is formed in 
the optically thick expanding stellar wind. In one par- 
ticularly remarkable display, the quintessential LBV S 
Doradus changed spectral type from B2e to F sometime 
between 1996 and 1999 (Masscy 2000). 

Our K band spectrum of WRAY 17-96 does not reveal 
an y significan t chan ges compared to the one presented 
in lEgan et"aLl ()2002f) , except in the highly variable He i 
2.058/Ltm line. To our knowledge our H band spectrum 
is the first one published of this source. Our K band 
spectrum of G026 .47-I-00.02 appears identical to that of 
iClark et all (|2003[ ). Our H band spectr um is of higher 
resolut ion and shows add itional lines not resolved in the 
one bv IClark et all ((200l . 

For G024. 73-1-00. 69, our spectra reveal dramatic dif- 
fer ences compa r ed to those from 2001/2002 presented 
by IClark et all (l200l . In our K band spectrum, the 
strength of the Na i doublet is significantly enhanced 
with respect to Br7, and the He i (2.112/zm) line has 
all but disappeared. In the H band, the changes are 
even more pronounced as the H I series has tran sitioned 
from emission to absorption. IClark et al.l (|2003[ ) derived 
a temperature of 12000 K from their 2001/2002 spectra 
of G024.73-f 00.69. However, the appearance of our spec- 
tra, in particular the absence of He i lines, now implies 
a significantly lower temperature. This spectral variabil- 
ity jJii_addition to the photometric variability discovered 
by IClark et al.l ()2003l ) , further solidifies the confirmation 
of G024.73-F00.69 as a bona fide LBV. 

The new LBV candidates identified in the previous 
section are strikingly similar in appearance to these 
known sources. Star #24 looks almost identical to #45 
(0026.47-^00.02). Star #14 resembles #29 (WRAY 17- 
96) with overall weaker features and star #17 is virtually 
a twin of #44 (G 024.73+00.69). Stars #46 and #23a 
stand out due to their pronounced [Fe ii] and Fe ii emis- 
sion, and bear a strong resemblance to the spectrum of 
WRA 751 in the K band (see Figure 8 of iMorris et al.l 
[1991) . 

Not surprisingly, our MIPS 24/im images of the shells 
(see Figure [T]) reveal significantly more detailed structure 
than the ISO images. The morphologies of the shells #29 
and #44 are remarkably similar, showing dense shells 
with fairly even surface bri ghtness. Both shel ls are also 
strongly detected at S/xm. IClark et al.l ()2003l) reported 
the possible detection of a much larger outer shell for 
#44. The reality of this outer shell is difficult to con- 
firm in our 8/im and 24/im images, given the complex 
structure in the general background emission evident at 
both wavelengths. It i s equally likely tha t the areas in- 
terpreted as "lobes" bv IClark et al.l ()2003[ ) are associated 
with unrelated embedded sources visible in the vicinity 
of that emission. In contrast to shells #29 and #44, shell 



#45 displays a more clumpy and filamentary structure. 
Some faint extended emission enhancements are seen at 
8/im coincident with the brighter areas at 24/xm, but the 
morphology of the extended emission close to the central 
star is quite different from that at 24/im. 

In addition to the known LBV / LBV candidates dis- 
cussed so far, shell #33 of our sample is associated with 
the LBV candidate HD 316285. We did not obtain any 
spectra of the cen tral source as part of this work, but 
iHillier et al.l (|1998| ) present optical and near-IR spectra 
of the star. Their H and K band spectra look almost 
identical to those of our central source #45. The star 
was known to be surrounded by a dust y shell based on 
the I RAS spectral energy distribution |McGregor et al.l 
Il988f ). The Spitzer 24/im image of HD 316285 is the first 
observation to actually resolve the shell. Its morphol- 
ogy, as well as Spitzer mid-IR spectroscopy of the shell 
is presented in .Morris (2008) . Note that there is some 
extended emission at 8/im as well, however, the morphol- 
ogy of that structure is very different compared to that at 
24/im. Overall, the characteristics of shell #33 are quite 
similar to those of shell #45. Amongst our newly iden- 
tified LBV candidates (#lla, 14, 17, 23a, 24, and 46), 
shells #14, #17, #24, and #46 most strongly resemble 
the shell morphologies of the known LBVs / LBV candi- 
dates, with shells #11 and #46 also exhibiting emission 
at 8/im. 

4.2. Optical Spectra 

Out of the nine central sources observed in the optical 
regime, two - stars #50 and #52 - are newly discov- 
ered WN stars (Figure [7]). We measured the equivalent 
width of N IV A4057, the N iii A4640 blend, and N V 
A4604, 4620 and followed t he classification scheme out- 
lined in iConti et al.l (|1990f ) . In addition, we compared 
our red spec t ra to those of the WR spectral catalog by 
iVreux et al.l (|1983f ). We determine a spectral type of 
WN6 for source #50. Note that this source was previ- 
ously classified as a planetary nebula, FN G029. 0-1-00. 4, 
its discovery dating back to 1966 (Abcll 1966). For star 
#52, we derive a subtype of WN7, confirming the clas- 
sification obtained earlier from our IR spectrum of this 
source. We also checked for the presence of hyd rogen fol- 
lowing the methodology of iSmith et al.l ()1996[ ). and ar- 
rive at a classification of WN6o for star #50 and WN7(h) 
for star #52. 

We utili zed the spectral atlases of [ Allen fc StromI 
(jl995) and iTorres-Dodgen fc Weaval ([19931) for an esti- 
mate of the spectral type for the remainder of the central 
sources. We also checked these spectral types through 
comparison with observations of spectral type standards 
obtained on the same nights as the targets. Source #59b 
is an early M type star, ~M2~M4 (the spectrum is not 
shown). The other six source spectra are displayed in 
Figure [5| In general, the blue spectra are quite noisy 
due to the highly reddened nature of these sources. Stars 
#39 and #44 were too faint to extract a signal at the 
shorter wavelengths, hence only red spectra could be ob- 
tained. Source #58, #60, and #61 are classified as B 
type stars (~B0-B5) based on the presence of He i (ab- 
sorptio n or emission) and absence of He ii A4541 and 
A4686 ([Walborn fc Fitzpatricklll990[ ). Source #60 has a 
pure H I, He i emission spectrum. Source #44 is the 
known LBV G024. 73-1-00. 69 discussed in the previous 
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section. Source #39 most likely also belongs to the group 
of B stars, but the red spectrum is quite noisy and it is 
difficult to gauge the reality of weak absorption features. 
Star #57 most closely matches an F/G type spectrum 
(~F5-G5). We are unable to determine the luminos- 
ity class for these sources based on our spectra. Note 
that the modulation in the continuum of the blue spec- 
tra (starting around 5300A) is an artifact of the data 
reduction process. 

5. INFRARED PHOTOMETRY 

2MASS and GLIMPSE photometry fo r each source has 
been asse mbled in Table [ 3 Fo llowing iMauerhan et al.l 
(|2009f) and lHadfield et al.lT2007D we constructed a color- 
color diagram (Figure |9]) to explore any correlations be- 
tween the known observed, and still unobserved sources 
of our sample. The comparison population of sources 
(black points) is composed of the photometrically most 
reliable point sources of a representative 1° x 1° "slice" 
of the Galactic plane from the GLIMPSE survey. The 
various types of sources we have identified spectroscopi- 
cally are indicated by different colors (see the caption of 
Figure [9] for details). 

Given the relatively small sample size for each group, 
one might expect limited insight from the distribution of 
points in color space, however, s ome general trend s seem 
to e merge. As already n oted in IMauerhan et ahl (|2OO90 
and lHadfield et all (|2007| ). the WR stars (red points) sep- 
arate quite cleanly from the main locus of field stars. The 
WR stars in our sample form a surprisingly well-defined 
and tight sequence, roughly parallel to the reddening vec- 
tor and offset to redder ifs— [8.0] colors from the bulk of 
the distribution of "normal" stars. The direction of the 
reddening vector at the locus of our WR stars is indi- 
cated by the dashe d line. Compared to t he analogous 
figure presented in IMauerhan et al.l (|2009l ). about half 
of our WR points have larger J — Ks values, indicating 
a more heavily reddened population. B stars and LBV 
candidates appear to populate the locus in between the 
normal stars and the WR track. Obviously, there is some 
overlap between the location of these different sources in 
the color-color diagram, since stars #24, #44, and #60 
lie on the sequence defined by the WRs in our sample. 
Note that there are also some outliers among the F-M 
group of sources that are on or close to the WR sequence 
(stars #23b, 26b, and 57). For the late type stars, this 
shift to redder colors might indicate a mid-IR excess due 
to cool circumstellar material. These objects could con- 
stitute a significant source of contamination for color- 
based searches for WR stars. 

Having identified the preferred loci for the various 
types of central stars, we now turn our attention to the 
as yet unclassified sources (green points). The reddest 
outlier at Xs— [8.0]> 4 is source #59a, the bluest at 
H — Ks K, 5.2 is source #27b. Their nature is unclear. 
Stars #41, #47, #48 and #54 are tightly clumped on 
the WR sequence at ifs— [8.0] « 2.1 and hence are the 
best candidates for WR stars waiting to be identified. 
Stars #2, #25, #34 and #55 are found among the main 
locus of field stars and most likely belong to the group 
of F-M type stars. The remainder of the sources of un- 
known type (with photometry in all bands) are sources 
#12, 15, 19, 27a, 30, 35, 38, 41, 47, 48, 49, 51, 53, 54 and 
56. We expect a mix of WR, LBV candidates, and F-M 



type stars among this group, but are unable to make any 
more finely differentiated predictions. 

6. SHELL MORPHOLOGY, SED, AND SIZE 

Morphologically, our 24/^m shells separate broadly into 
two categories: fainter, more elliptical shells with some- 
what irregular rims, and strikingly circular nebulae with 
sharply defined edges. Shells associated with late type 
stars (our F-M types) overwhelmingly belong to the first 
category (with the exception of shells #8 and #57), while 
the early type stars are usually exhibiting the most cir- 
cular shells (exceptions are #16, 18, 25, and 32). No 
distinguishing characteristics are evident in the appear- 
ance of the shells between WR stars, LBV candidates 
and B stars. The morphology of some of the shells sur- 
rounding known LBV stars has already been discussed in 
Section 14.1.31 Based on the general morphological dis- 
tinction at 24/im between late and early type stars, we 
conclude that the hot star is likely producing the shell 
in the case of the double central sources containing both 
a late and early type star (#23, #28). We also note 
that the morphology of shell #52 provides evidence for 
multiple, discrete mass loss events. Several concentric 
rings can be discerned even in the presence of significant 
structure in the overall background emission. 

As noted in Section 1, many of the shells detected at 
24/im do not exhibit equivalent emission at 8/im: 37 out 
of the 61 shells listed in Table [1] (60%) would not have 
been identified as shell sources at 8/im. Considering each 
of the spectral type subcategories, both of the known 
PNe, five out of the nine late type stars, none of the B 
type stars, two out of the nine Oe/WN and five out of 
the 10 Be/B[e]/LBV type stars show 8/im shell emission. 
The mid-IR spectroscopic observations of some LBV and 
WR shells imply that the non-detection of 8/im emission 
could simply be explained by a reduction in flux density 
at this waveleng th relative to that at 24 /im. The spec- 
tra pr esented by iVoors et al.l ()2000[ ) and iBarniske et ahl 
(|2008f) show a steep increase in emission towards longer 
wavelength starting at around 10/im. At the same time, 
we cannot exclude the possibility that emission only ob- 
served at 24/im might be caused by a strong line feature 
in the spectrum, s uch as the [OIV] li ne producing the 
shell discovered bv lMorris et al.l ()2006[ ). 

Based on the general morphological properties of the 
as yet unclassified sources, we speculate that the shells 
#2, 12, 25, 34, 59a are formed by late type stars. Recall 
that we predicted sources #2, 25, 34, and 55 to likely be 
stars of type F-M from their location in the color-color 
diagram, supporting this conclusion. We predict that 
shells #7, 27, 30, 35, 37, 38, 41, 47, 48, 49, 51, 54, and 56 
are associated with early type central stars. Shells #15, 
19, 40, 42, 53, 55 remain ambiguous, partially because 
none of these have a prominent central source at 24/im, 
unlike all of the other shells. 

The observed radii of the 24/im shells listed in Table [T] 
range from 0.14-2.4' with an average value of ~ 0.7'. The 
physical sizes for the shells surrounding late type stars 
have been derived in Section 14.1.1! and are listed in Ta- 
ble [21 We now also estimate the dimensions of the shells 
around WR type sources among our sample, adopting 
the intrinsic colors and absolute Ks ban d magnitudes for 
the various WR subtypes presented by iCrowther et al.l 
((2006|). Utilizing the extinction ratios for the JHKs 
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bands from llndebetouw et all (|2005f) , we first calculate 
an average extinction Aks based on the observed 2MASS 
J — Ks and H — Ks colors (provided in Table [3]) . This 
extinction is combined with the observed Ks magnitude 
and the assumed M^s for the appropriate WR subtype, 
and we thus arrive at the distances and physical shell 
sizes listed in Tabled Note that in addition to our newly 
identified WR stars #10, #l lb, #50 and #52, we also 
list t he known WR stars #9 (|Mauerhan et al.ll2009[ ) and 
#16 (|Shara et al.ll2009t ). The shehs around these two 
sources have not been noticed previously, so the physi- 
cal size estimate for the shells are new calculations, while 
the distances and Aks values for these stars were already 
presented in their respective discovery references. The 
WC9 star #llb is one of two possible central sources for 
shell #11 (the fitted center of the shell lies in between 
the positions of the two stars; see Section 8). We clas- 
sify the other star, #lla, as belonging to the group of 
Be/B[e]/LBV type stars, hence either star could be pro- 
ducing the observed shell, although circumstellar nebu- 
lae are more common to B[e]/LBV stars than the more 
highly evolved WC stars. The close association of two 
such rare sources makes it likely that they are located at 
the same distance and could indicate the presence of a 
small cluster of stars. An investigation of the surround- 
ing sources might be worthwhile. In any case, a distance 
estimate for one of the sources will deter mine the physica l 
size of the shell. Unfortunately, Cro wther et al.l ()2006f ) 
do not provide a Mks for WC9 stars, only intrinsic IR 
colors. We have therefore adopted the M^s magnitude 
of star E from their Table 9 as a rough estimate. 

The extinction, distance and size scales listed in Ta- 
ble H] seem reasonable, except in the case of star #50, 
where a large distance (d=16 kpc) is paired with a rela- 
tively low extinction value {Ay — 5.7). This is one of the 
few objects that exhibits shell emission at optical wave- 
lengths. As already mentioned, the shell was classified 
as a planetary nebula, but we were unable to find any 
previous classification of the central star (see Section |4?2] 
and 8). 

Comparing the radii of our WR shells to those reported 
in the literature, our shells generally fit well within the 
observed range, possibly somewhat skewed to the smaller 
end of the distribution . An optical study of eight WR 
ring nebulae by Grucn dl et al.l (|2000D finds radii between 
2.1-9.0 pc, and iChul (|1991| ) quotes a few to 30 pc as the 
typical diameters of WR shells. [Marston (1996) surveyed 
156 WR stars with IRAS and finds 49 "probable or sus- 
pected shells" ranging in size between 5-150 pc in radius. 
Due to the low spatial resolution of IRAS, these detec- 
tions are clearly biased towards the largest shells. 

We are unable to determine the physical sizes 
for the remainder of our shells, mostly comprising 
B/Be/B[e]/LBV type stars. It is problematic to derive 
accurate distances for field B supergiants because they 
lack luminosity diagnostics in their spectra. Detailed 
modeling of the stellar IR spectra as well as the shells 
to constrain the luminosity and reddening of the stars 
and determine shell parameters such as mass and dust 
composition, are beyond the scope of this paper and will 
be addressed in future work. iClark et al.l ()2003[ ) summa- 
rize the best estimates for the radii of LBV and candidate 
LBV shell nebulae, which range from 0.1-3.6 pc, most of 
them smaller than ~ 1 pc. 



7. SUMMARY AND DISCUSSION 

We have discovered a large number of prominent, 
highly symmetric 24/im shells surrounding bright central 
sources. We were able to determine the nature of the ap- 
parent central sources for 45 of the 62 shells, 10 based on 
literature searches and 35 from dedicated spectroscopic 
follow-up. Our infrared and optical spectroscopy has re- 
vealed a mix of early and late type stars. However, the 
prevalence of rare, massive stars stands out. We find 
six bona fide WR stars in our sample (four of which are 
new identifications), as well as four stars which we ten- 
tatively classify as Oe/WN. Most strikingly, our sample 
of shells contains four known LBV type stars (one con- 
firmed, three candidates), and we classify an additional 
six sources as LBV candidates based on the similarities 
between the spectra. We realize that the classification of 
a particular source as an LBV involves a host of crite- 
ria, including photometric and spectroscopic variability, 
that still need to be investigated for these new sources. 
However, given the small numbers of known LBVs, the 
presence of a shell combined with the spectroscopic char- 
acteristics make this relatively large sample of new can- 
didates an exciting discovery. We also identify five new 
Be stars in our sample. The distinction between these 
and what we characterize as LBV candidates is some- 
what arbitrary and simply conveys less peculiar spectral 
properties. It is entirely possible that these stars harbor 
additional LBV candidates, especially considering LBV 
spectral variability and the fact that the LBVs AG Car 
and HR Car have been classified as Be stars. 

The detailed evolutionary sequence for the most mas- 
sive stars is still a matter of intense debate. The standard 
model of massi ve star evolution (see, e.g., the review by 
lCrowthed2007| ) expounds that O stars with initial masses 
of 25-40Mq will undergo an LBV and RSG phase on 
their way to becoming a WR star. More massive stars 
on the other hand, are thought to avoid the RSG phase. 
Smith & Owocki (2006) contend that substantial mass 
loss during an LBV stage is critical to the formation 
of WR stars. Analysis of the nebulae surrounding the 
known LBVs support the view that the n ebulae formed 
durin g a BSC and not d uring a RSG phase ()Lamers et al.l 
I2001D . iSmith fc Contil (2008) strongly argue in favor of 
the notion that H-rich WN type stars precede the LBV 
stage, while H-poor WN stars represent a post-LBV evo- 
lutionary phase. In this context we point out that our 
Oe/WN classified sources all show strong H emission, and 
hence could be interpreted as pre-LBV stars. In fact, we 
speculate that a large number of the central stars associ- 
ated with our most well-defined, circular shells represent 
a narrow evolutionary phase, closely related or identical 
to the LBV phase. In depth analysis of the shell emission 
is needed to further investigate this hypothesis. Finally, 
we predict that our upcoming observations of the remain- 
der of our sample will reveal additional members of this 
pre-LBV/LBV group among the as yet unclassified shell 
central stars. 

8. APPENDIX: NOTES ON INDIVIDUAL OBJECTS 

1: IRAS 11419-6228 is coincident with the central source 
position. 

2: M O 1-184, class ifi ed as a OB- star by 
iMuzzio fc Orsattil (|1977l) . is offset by 35.18" 
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from our central source position and has a 24/im 
counterpart at the edge of the shell, coincident 
with an emission enhancement in the shell. Hence 
MO 1-184 is distinct from the central source 
associated with shell #2. Note that the 24/im 
central source is offset from the IRAC 3.6 and 
8.0/im source. In the 2MASS images, several 
blended sources are visible at the central source 
position. The 2MASS source listed in Table [1] 
corresponds to the source that coincides with the 
brightest IRAC source within 1" of the center 
position. 

4: IRAS 14173-6124 is offset by 21" away from the cen- 
tral source, but is likely associated with the central 
source. At 2MASS wavelengths there is a group 
(cluster?) of bright stars that coincides with the 
southeastern border of the shell. 

5: IRAS 15039-5806 is 64" away and clearly associated 
with a different bright 24/im point source. In 
2MASS, several sources are blended at the central 
source position. 

6: The fitted geometric center of the ellipse falls halfway 
between the two bright 24/im sources. 

8: IRAS 15293-5602 is offset by 32" from the central 
source and clearly associated with a bright 24/im 
source on the east rim of the shell. There are four 
bright 24/im sources within the shell, source #8a 
is at the geometric center, but we also obtained 
a spectrum of the source immediately to the west 
(#8b). 

9: This WR star was identified based on 2MASS 
and GLIMPSE photo metric selection criteria by 
iMauerhan etall pOOl ). The presence of the shell, 
however, has not been noticed previously. 

10: IRAS 15421-5323. 

11: 2MASS J15491137-5508516 IRAS 15452-5459), 
260" distant, is classified as a post AGB star 
(proto-PN) in SIMBAD. It is associated with an 
extremely bright 24/tm point source i n the MIPS- 
GAL images, but not with our shell. ISahai et al.l 
()2007D present HST imaging of the post- AGB star 
which clearly shows a double lobed structure. How- 
ever at 24/im there is no large scale extended emis- 
sion associated with that source, instead only our 
shell stands out. The fitted geometric center of the 
circular shell falls halfway between the two 24/im 
sources listed in Table [TJ 

12: IRAS 15517-5334 is 40" offset from our shell central 
source, but lies within the shell. The association is 
uncertain. 

14: IRAS 16254-4739 (15"offset) clearly associated with 
our shell. The only refere nces for this object 
are Ivan der Walt et al.l (|1996l ) and lBronfman et all 
(199^ who search for 6.7 GHz methanol maser 
emission associated with IRAS sources selected to 
meet criteria for being ultracompact H ii regions. 
No maser emission was detected for IRAS 16254- 
4739. 



15: An association with IRAS 16278-4808 (28" away) 
is possible, but somewhat ambiguous. C* 2338, a 
carbon star with a position offset by 122" is clearly 
associated with another 24/tm source. 

16: IShara et al.l (l2009t ) list a WN5b star (Shara 1093.53) 
near the center of this shell, which is difficult to pin- 
point due to its diffuse outer edge. They did not 
detect the shell itself, and there are three bright 
sources within about 8" of the central region. This 
is one case where selecting the source coincident 
with the brightest central 24/im source would have 
identified the wrong star, as the WN5 star is not 
associated with any enhanced 24/im emission. We 
note that IGR J16320-4751, an embedded high 
mass X-ray binary, is also found in the vicinity 
(157" offset). 

18: IRAS 16396-4555, EM* VRMF 115, Hen 2-179, 
WRAY 16-232, SS73 63, has been classified as any- 
thing from PN to "not a PN" , highly obscured Be 
star to M supergiant, or peculiar emission line star. 
A good summary of the literature on this source, 
as well as an optic al spectrum, are prese nted in 
IPereira et all (|2003D . IPereira et all (|2003l ) classi- 
fies the source as a reddened Be star. Note that 
the shell has a highly clumped morphology. 

19: IRAS 16426-4504 (60" offset) is clearly associated 
with another 24/im source. IGR J16465-4507 is 
located nearby (215"), but not coincident with our 
shell. The IGR source consists of a compact object 
and a blue supergiant companion of spectral type 
09.5Ia (Nespoli et al. 2008). 

21: The planetary nebula IC 4637 is discussed as a pos- 
sible containing a physical double central source 
with a K4 com p anion at a separation of 2.42" by 
iCiardullo et"aLl (|1999[) . 2MASS clearly shows two 
sources, but only one appears in the 2MASS Point 
Source Catalog. The GLI MPSE morphology of the 
PN sh ell is described in iPhillips fc Ramos-LariosI 
(l2008l) . 

22: IRAS 17039-3952 is offset by 30" located halfway 
between our shell and another embedded source. 
The association is ambiguous. 

23: IRAS 17050-3921. 

26: HD 317949, classified as a KO star in SIMBAD, is off- 
set by 36" and definitely associated with another 
source at 24/i,m at the edge of the shell. The ge- 
ometric center of the shell is located between the 
two 24/im sources. 

27: =i^27a is the source closest to the geometric center 
of the shell. At 2MASS wavelengths, the central 
two sources appear to be blended with additional 
sources. 

28: The ra dio source GPA 357.18-0.14 (|Langston et al.l 
I2OOO0 is found at a separation of 24". It is not 
associated with our central source, but located on 
the SW rim of the shell. Source #28a is closest to 
the geometric center of the shell. 
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29: The LBV candidate WRAY 17-96 was originally 
thought to be a PN. lEgan et alj |2002) describe 
the detection of a ring nebula with MSX and show 
a K band spectrum of the central source. 

31: Planetary Neb ul a PN G003.5+02.7, discovered by 
iBoumis e t al.l (I2003D and in dependently by the 
MASH survey ()Parker et al.l [2006. ). Morphologi- 
cally, the 24/im shell strongly resembles the Ha 
shell from the MASH survey. Our identification of 
the central source is uncertain due to a slight offset 
of the selected star from the geometric center of the 
shell. 

32: IRAS 17408-3027. IRaharto et all (IT981 list [RHI84] 
10-469 at the position of our central source and give 
a M2: type. No additional information is provided. 

33: H P 316 285 is a known LBV candidate. IHillier et afl 
(!1998D present optical and near-IR spectra of the 
central source. The star is known to be surrounded 
by a dusty sh ell based on the I RAS spectral energy 
distribution ()McGregor et al.l 11988). The MIPS 
24/im image is the first observation to resolve the 
shell. Note that there is some extended emission 
at 8/im as well, however, the morphology of that 
structure is very different from that at 24/im. 

36: At 2MASS wavelengths, the central source is clearly 
a blend of two sources. However, only one source 
is listed in the 2MASS catalog. The location of 
the source is within the coverage area of the MAG- 
PIS survey (jHelfand et al.ll2006[) . however, no ra- 
dio emission is detected from either the shell or the 
central source. 

41: The X-ray source XGPS-I J182833-102652 is located 
118" away, but clearly associated with a different 
24/im source and nebulosity. At 2MASS wave- 
lengths, the central source is blended with at least 
one other star. The location of the source is within 
the coverage area of the MAGPIS survey, however, 
no radio emission is detected from either the shell 
or the central source. 

44: G 024.73-H00 . 69 (= V481 Set) is a known LBV 
(jClark et al.ll2005D wit h a known shell. A K band 
spectrum is shown in iClark et"al1 (|2003[ ) . No ra- 
dio emission is detected from either the shell or the 
central source with the MAGPIS survey. 

45: G 26.470-FO 021 is a known LBV candidate 
(jClark et al.l l2005l with a known shell. The MAG- 
PIS 20cm image shows bright emission coincident 
with the central source, as well as some nebular 
structure. 

47 and 48: Source ^^47 and #48 form a spectacular 
complex of two neighboring shells at 24/im. We 
treat each shell component as a separate entity, al- 
though the combined structure is also reminiscent 
of a double lobed nebula. However, there is no 
obvious central source at the "waist" of this po- 
tential bipolar emission. It is unclear whether the 
two shells are interacting with one another or sim- 
ply are a chance superposition of two independent 



shells. Due to the large size of each shell, the identi- 
fication of the central sources is uncertain. At radio 
wavelengths, the outline of shell #48 is faintly visi- 
ble in the MAGPIS 20cm image, as well as a bright 
point source near the southern edge of the shell. 

49: Both the shell and the central source are visible in 
the MAGPIS 20cm image. 

50: This source is identified as planetary nebula PN 
G029.0-H 00.4 i n SIM BAD. Its discovery dates 
back tolA belll C 1 9661) wh ere it is listed as A48. 
ICondon fc KaplanI (|l998l ) include the object in 
their catalog of PN in the NRAO VLA sky survey, 
but note it as "Confused by complex background on 
the NVSS image." We were unable to find any pre- 
vious classification of t he central star, e.g., no cen- 
tral star type is listed in lZuckerman &: Alleii 0986). 
The shell, but not the central star, is a strong radio 
source and clearly defined in the MAGPIS image. 

51: SIMBAD lists [BSM2002] 18431-0312 4, a cloud of 
unknown nature at a distance of 53". This is also 
IRAS 18431-0312 which in the 24/im MIPSGAL 
image is clearly associated with a different object, 
distinct from our shell. No radio emission is de- 
tected from either the shell or the central source 
with the MAGPIS survey. 

56: IRAS 18588-1-0350 is offset by 43" the association is 
uncertain. No radio emission is detected from ei- 
ther the shell or the central source with the MAG- 
PIS survey. 

59: Source #59a is closest to the geometric center of the 
shell. No radio emission is detected from either 
the shell or the central source with the MAGPIS 
survey. 

60: This source is listed as IRAS 19425-^2411 and 
HBHA 2203-01 in SIMBAD The only reference 
([Kohoutek fc Wehmeveii [19991) simply includes the 
central source in a catalog of Ha emission line stars. 
No additional information is provided. 

62: The central star of this known PN, NGC 6842, is a 
blend of two sources at 2MASS wavelengths. 



This work is based in part on archival data obtained 
with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute 
of Technology under a contract with NASA. Support 
for this work was provided by an award issued by 
JPL/Caltech. Based on observations obtained at the 
Hale Telescope, Palomar Observatory, as a part a con- 
tinuing collaboration between the California Institute of 
Technology, NASA/JPL, and Cornell University. This 
publication makes use of data products from the 2 Micron 
All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analy- 
sis Center /California Institute of Technology, funded by 
the National Aeronautics and Space Administration and 
the National Science Foundation. It also utilized NASA's 
Astrophysics Data System Abstract Service and the SIM- 
BAD database operated by CDS, Strasbourg, France. 
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TABLE 1 

Shell Data 



Num 2MASSJ Designation radius central src 


shell at 


Comments^ 


(arcmin) 24/8/3.6/1.1 


8^m? 





1 


1 1441oUo-DZ4DZlU 


U.Oo X U.vSO 


y/y/y/y 


no 


o 


1oU1z1oo-DZ441o1 


n ri/i 1 ni 
U.y4 X l.Ul 


y/y/y/y 


yes 


Q 
O 


lolU4oo4-Dol ( 40 i 


n A^ 

U.40 


y/y/y/y 


no 


4 


14zilUUZ-Div5oiyy 


1 no ,r o Qo 
l.Uo X Z.oo 


y/y/y/y 


yes 


5 


loU/ODzl-oolozoU 


n c;n -.r n /I o 
U.OU X U.4Z 


y/y/y/y 


no 


ua 


iouyzoyy-oo4 / 4ou 


n An n on 
U.4y X U.oU 


y/y/y/y 


yes 


t)D 


iouyzyoi-0o4 / ooi 




y/y/y/y 




( 


lOZOOyou-O / U44Uv5 


n 1 7 
U.l ( 


y/y/y/y 


yes 


fin 

oa 


lOvSoUouy-ooizzuu 


U.DD 


y/y/y/y 


no 


QK 
OD 


lOvSoUDDo-ODlZZiy 




y/y/y/y 




y 


10v50ZdOZ-OdU41Zv5 


1 n7 
l.U / 


y/y/y/y 


no 


1 n 


10400yi4-0ov5Zv5Z0 


n 7n 
U. /U 


y/y/y/y 


no 


1 1 «^ 
lia 


104o4ZU /-OOU / 4ZZ 


1 c;7 
1.0 / 


y/y/y/y 


no 


1 1 K 


104o4ZlU-OoU / 04Z 






12 


lOOOO / OO-0o4o4Uz 


1 O /I -.^ O /I o 

i.z4 X Z.4Z 


y/y/y/y 


yes 


1 

lo 


lOOolo { o-OZO / Olo 


U.oO 


y/y/y/y 


no 


1 A 


IDZyUoi /-4/4uZ04 


U.OO 


y/y/y/y 


yes 


10 


lUOlO / Ol-4ol400o 


U.ol X U.Zo 


i^/y/y/y 


no 


16 


iDoZizyo-4/oUoOo 


1.02 


n/y/y/y 


no 


17 


10v5o4Z i 0-4oouzU ( 


0.30 


y/y/y/y 


no 


1 Q 

lo 


lu4oluoD-4uUU4Z4 


u.yo 


y/y/y/y 


no 


1 o 

ly 


lu4ul /o4-40Uo4/O 


n 1 a 
U. lo 


y/y/y/y 


no 




iu4yo/ (U-4ooooyz 


1 /I a 

1 .4o 


y/y/y/y 


no 


zi 


1 / U01U4o-4U0oU / 1 


n oc; 

U.ZO 


ii/n/y/y 


yes 


zz 


1 / U / Zv>v)o-oyOuOU4 


n /iQ n 71 
U.4o X \J. 1 1 


y/y/y/y 


yes 


zoa 


1 /Uozyl^>-oyzou/f^ 


n 7n n 7*3 
U. i y X U. r v> 


y/y/y/y 


no 


ZoD 


1 / uozyou-oyzDioo 




y/y/y/y 




0/1 


1 / iiuuy4-oy4Di / 4 


n Ad 
U.4u 


y/y/y/y 


no 


oc; 
ZO 


1 / oOZy4Z-oZ4oOuZ 


U.OO X U.4/ 


y/y/y/y 


no 


26a 


1 (o(\j6 I l-ol4/4o ( 


1 nn -.^ n no 

i.uy X u.yz 


y/y/y/y 


yes 


26b 


I/O / Uoy4-ol4 / o4y 




y/y/y/y 




2 ra 


1 / v5/4/o4-olv5/oov5 


0.51 


y/y/y/y 


yes 


27b 


1 / o /4/oU-olo (ofu 




n/y/y/y 




OQo 

Zoa 


1 7'?01 SOO QTO/IOQO 

1 /oyioyy-oiz4Zoy 


n KK 
U.OO 


y/y/y/y 


no 


28b 


i / oyiyio-v5iz414z 




y/y/y/y 




29 


17/11 Qc;/iQ onrii^ q on 
1 / 4ioo4o-oUUuooy 


1.0 


y/y/y/y 


yes 


oU 


i7/ioi/im onc^cQcn 
1 / 4Z14Ul-ZyOOoOU 


n OB 
U.Zo 


y/y/y/y 


yes 


Q 1 

ol 


i7/i'j'3nnQ 0/iQic;oc; 
1 / 4^>oyUo-Z4^510Z0 


n o/i 
U.Z4 


i^/y/y/y 


yes 


QO 


1 / 4o0yol-OUZoOO4 


1 r\A V n 7fi 
1 . U4 xU.ru 


y/y/y/y 


no 


QQ 
OO 


17/iQi/in'j OQnnc;Qi 
1 / 4ol4U^>-ZoUU0ol 


1 c;/l 1 n7 
1.04 X i.y / 


y/y/y/y 


yes 


34 


1 / 4yoOv5y-zo4yoUz 


n £;o -.^ n 7n 
U.bo X U. r y 


y/y/y/y 


yes 


35 


1 onn c; 7^? o o /i /i ^; 7 
loUUo / DZ-24oo4d ( 


0.17 


n/y/y/y 


yes 


OD 


1 snoooQ'i oo'^anno 
Lo\JZZZoo-zzoo\jyjz 


n /io 

U.4Z 


y/y/y/y 


yes 


O i 


loUoOOO / -ZZOuUUU 


n Ol 
U.Zl 


y/Vy/y 


yes 


OO 


loU / Uolo-zulDluo 


U.OU 


y/y/y/y 


no 


oy 


isi'i'iioi 1 sc;i?/iQi 
Ioiooizi-io0u4oi 


n oa 
U.Zo 


n/y/y/y 


no 




ifiOQQ/ion 1 ciOE^n7fi 

loZoo4ZU- lOZOU ( D 


n 07 

U.Z i 


n/y/y/y 


no 


/1 1 
41 


1oZo410D-1UZ ( UOD 


n 77 
U. / ( 


y/y/y/y 


no 


/I o 
4Z 


1 cioi 7*50 noi ^^1 Qo 
lovjzi ( vSy-uyiDioy 


n 1 /i 

U. 14 


n/n/y/y 


no 


4v5 


Iovjooy04-Uou ( Uo4 


U.OO 


y/y/y/y 


no 


/I /I 
44 


loooOOZo-UoOovSoD 


n 7/1 
U. ( 4 


y/y/y/y 


yes 


40 


1 0*50 '500/1 nc;/i /I on/i 

lOv5yv5ZZ4-U044ZU4 


1 '30 V 1 n/1 
i.vSy X 1.U4 


y/y/y/y 


yes 


4D 


lo410yDD-U0l04Uy 


n A'i 

U.4o 


y/y/y/y 


yes 


A 1 
4 ( 


IQ/IOrii^QH riQ/IQOO/l 

1o4ZUDoU-Uo4oZZ4 


1 fi7 
l.D ( 


y/y/y/y 


no 


A Q 
4o 


iQ/ionfi07 riQf^inoQ 
io4ZUoz i-uoOiuzy 


1 If; 

1.10 


y/y/y/y 


no 


4y 


1Q/1000/17 nc;n/iQnn 
lo4ZZZ4 /-U0U4oUU 


n /1 7 
U.4 / 


y/y/y/y 


no 


50 


18424692-0313172 


0.29 


n /v /v /v 


yes 


51 


18455593-0308297 


0.40 


y/y/y/y 


yes 


52 


18492733-0104207 


1.8 


y/y/y/y 


no 


53 


18503980-^0004453 


0.33 


n/y/y/y 


no 


54 


18510295-0058242 


1.07 


y/y/y/y 


no 


55 


18530582-^0011358 


0.24 


n/y/y/y 


no 


56 


19011669-^0355108 


0.52 


y/y/y/y 


yes 


57 


19042098-^0600001 


0.74 


y/y/y/y 


no 


58 


19325284+1742303 


0.78 


y/y/y/y 


no 


59a 


19385510-^2127550 


0.89 X 0.49 


y/y/y/y 


yes 


59b 


19385569-^2127483 








60 


19443759-^2419058 


0.53 


y/y/y/y 


no 


61 


19444295+2311337 


0.46 


y/y/y/y 


no 


62a 


19550232+2917178 


0.42 


n///y 




62b 


19550249+2917198 




y///y 





IRAS 11419-6228; Oe/WN (IR), this work 
Oe/WN (IR), this work 

IRAS 14173-6124; M6 III / K5 I (IR), this work 
MO III / KO I (IR), this work 
M6 III / K5 I (IR), this work 
M5 III / K4 I (IR), this work 

Kl III / G7 I (IR),this work 

M5 III / K4 I (IR),this work 

WN7 (Maue rhan et ah 2009) 

IRAS 15421-5323; WNghTTR), this work 

Be/B[c]/LBV (IR), this work 

WC9 (IR), this work 

IRAS 15517-5334: 

Oe/WN (IR), this work 

IRAS 16254-4739; Be/B[e]/LBV (IR), this work 
IRAS 16278-4808: 
WN5b (Shara ot al. 2009) 
Be/B[e]/LBV (IR^" this work 

IRAS 16396-4555, SS 73 63; Be (Pereira et al. 2003) 
Be (IR), this work 

PN G345. 4-1-00.1, IC 4637, Hen 2-193, WRAY 15-1607 
Ml I (IR), this work 

IRAS 17050-3921; Be/B[e]/LBV (IR), this work 
M2 III / K2 I (IR), this work 
Be/B[e]/LBV (IR), this work 

MO III / KO I (IR), this work 
G9 III (IR), this work 



OB (IR), this work 

M6 III / K5 I (IR), this work 

WRAY 17-96, Hen 3-1453, LBV candidate, (IR) 

PN G003. 5-1-02. 7 

star M2:; Be (IR), this work 

HD 316285, Hen 3-1482, WRAY 15-1777, LBV candidate 

IRAS 17578-2433 
Oe/WN (IR), this work 



Be: (OPT), this work 



G7 I (IR), this work 

G 024.73-f00.69, V481 Set, LBV, (IR, OPT) 

G 026. 47-f 00.02, 2MASS J18393224-0544204, LBV candidate, (IR) 
Be/B[e]/LBV (IR), this work 



IRAS 18397-0507 

PN G029.0+00.4; WN6 (OPT), this work 

WN7 (IR, OPT), this work 
IRAS 18481+0001 



IRAS 18588+0350: 
F/G (OPT), this work 
B0e-B5e (OPT), this work 

M2-M4 (OPT), this work 

IRAS 19425+2411, HBHA 2203-01; Be (OPT, this work) 

B0-B5 (OPT), this work 

PN G065.9+00.5, NGC 6842, Hen 2-451 



^ For detailed notes on each object see Section 8 
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TABLE 2 
Late Type Sources 



Source 


Type 


Av 


d 


size 






(mag) 


(kpc) 


(pc) 


4 


M6 111 


29.8 


1.1 


0.34 X 0.75 




K5 1 


32. -5 


2.6 


0.80 X 1.77 


5 


MO 111 


7.8 


1.2 


U.17 X U.i4 




KO 1 


10.4 


6.1 


0.89 X 0.75 


6a 


M6 111 


12.0 


3.2 


0.46 X 0.75 




K5 I 


14.7 


7.6 


1.08 X 1.76 


6b 


M5 III 


13.8 


2.5 


0.36 X 0.58 




K4 I 


16.7 


7.6 


1.U9 X 1.77 


!5a 


l^n TTT 
lv± 111 


T A A 
±4.4 


U.O 


nil 
U.ii 






15. o 


9.6 


1.85 


8b 


M5 III 


12.2 


4.1 


0.78 




K4 I 


15.0 


12.5 


2.39 


22 


Ml I 


8.6 


4.3 


0.53 X 0.88 


23bl 


M2 111 


13.1 


2.7 


0.63 X 0.58 




K2 1 


15.7 


11.6 


2.67 X 2.47 


26a 


MO III 


10.2 


2.1 


0.67 X 0.56 




KO I 


12.8 


11.0 


3.50 X 2.96 


26b 


G9 III 


9.6 


1.4 


0.44 X 0.37 




G5 1 


10.4 


23.7 


7.51 X 6.34 


28bl 


M6 111 


18.2 


5.3 


0.86 




K5 1 


20.9 


12.6 


2.01 


43 


Kl III 


24.8 


0.3 


0.44 




G7 I 


25.6 


4.6 


0.74 


Probably not producing the shell, see 


text for details. 
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TABLE 3 
Shell Central Source Photometry 



Num 


OTV IT A CO T 

2MAbb J 


tW T A CO TT 


oil (T A CO 


3.6/im^ 


Ac: 1 


CO 1 

5.8^m'^ 


8.0/im^ 


1 


9.999(23) 


8.720(31) 


7.775(23) 


6.982(26) 




6.175(37) 


5.813(30) 


2 


15.090(66) 


12.459(75) 


11.221(36) 


10.281(31) 


10.250(42) 


10.001(36) 


10.085(33) 


3 


12.271(23) 


11.004(20) 


10.186(19) 


9.516(31) 


9.065(41) 


8.824(28) 


8.450(22) 


4 


12.719(48) 


8.618(55) 


6.448(15) 


5.151(112) 




4.288(29) 


4.169(29) 


5 


9.282(20) 


7.815(21) 


6.959(13) 


6. 600(61 ) 


6 1AA(57) 


5.900(37) 


5.882(33) 


6a 


10.061(19) 


7.905(39) 


6.776(7) 


6. 697(74 ) 




5.731(33) 


5.594(26) 


6b 


10.870(23) 


8.520(41) 


7.311(9) 


6.617(46) 


6.482(46) 


6.142(35) 


6.020(27) 


7 


15.759(94) 


13.628(77) 


12.390(57) 


11.929(309) 


11.851(194) 


11.280(97) 




8a 


11.678(28) 


9.688(20) 


8.573(19) 


7.572(36) 


7.274(47) 


7.043(32) 


6.890(27) 


8b 


11.472(21) 


9.218(20) 


8.189(29) 


7.510(47) 


7.592(46) 


7.330(34) 


7.228(23) 


9 


13.841(35) 


12.386(39) 


11.457(27) 


10.583(49) 


10.120(59) 


9.843(47) 


9.490(31) 


10 


14.576(34) 


12.217(36) 


10.874(21) 


9.447(64) 


9.023(55) 


8.620(48) 


8.399(29) 


11a 


7.929(23) 


6.725(45) 


5.977(13) 




5.142(68) 


4.943(26) 


4.783(27) 


lib 


11.185(25) 


9.708(28) 


8.490(21) 


7.203(145) 


6.632(42) 


6.239(37) 


6.024(28) 


12 


12.380(25) 


8.637(17) 


6.603(17) 


5. 215(182) 




3.923(25) 


3. 969(33) 


13 


13.774(25) 


11.369(29) 


9.984(27) 


8.652(44) 


8.185(39) 


7.854(37) 


7.650(34) 


14 


8.668(27) 


7.141(43) 


6.188(15) 


6.749(106) 


5. 626(63 ) 


5.076(31) 


4.735(23) 


15 


>16.567 


12.595(25) 


10.281(18) 


8.173(46) 


7.462(40) 


7.107(30) 


6.854(30) 


16 


14.998(52) 


12.697(33) 


11.340(26) 


10.049(49) 


9.514(44) 


9.196(88) 


8.988(46) 


17 


11.711(18) 


7.889(39) 


5.820(15) 


A 915fll ) 


4.519(34) 


4.046(17) 


4.054(15) 


18 


6.258(17) 


5.077(45) 


4.205(33) 










19 


15.380(102) 


11.867(30) 


9.847(21) 


8.091(36) 


7.562(41) 


7.252(36) 


7.148(32) 


20 


7.242(17) 


6.086(33) 


5.422(17) 




5.022(55) 


4.518(32) 


4.336(19) 


21 


11.364(61) 


11.343(91) 


10.964(63) 










22 


7.009(9) 


5.325(13) 


4.635(13) 






4.096(28) 


4.109(23) 


23a 


10.492(21) 


8.219(51) 


7.111(19) 


6.463(87) 


6. 1 64 (55 ) 


5.714(26) 


5.545(29) 


23b 


12.081(21) 


10.032(24) 


8.790(17) 


7.686(47) 


7.151(53) 


6.788(32) 


6.442(25) 


24 


9.728(21) 


8.041(33) 


6.963(17) 


6.673(154) 


6.075(88) 


5.239(37) 


4.884(23) 


25 


10.749(19) 


9.101(20) 


8.254(17) 


7.633(36) 


7.443(29) 


7.295(29) 


7.336(19) 


26a 


11.260(23) 


9.342(23) 


8.527(21) 


7.949(31) 


8.093(29) 


7.885(30) 


7.872(19) 


26b 


12.556(23) 


11.014(26) 


10.339(32) 


9.365(49) 


9.038(90) 


8.711(34) 


8.536(21) 


27a 


11.771(34) 


9.248(38) 


7.788(9) 


6.838(42) 


6.191(45) 


5.957(28) 


5.834(22) 


27b 


13.701(6) 


10.588(41) 


8.464(49) 


8.095(99) 


8.316(89) 


7.747(42) 


7.921(38) 


28a 


8.440(19) 


7.586(29) 


7.078(11) 


6.864(33) 


6.657(44) 


6.548(27) 


6.573(23) 


28b 


>12.898 


10.078(36) 


8.578(31) 


7.464(28) 


7.502(29) 


7.080(28) 


7.029(18) 


29 


6.707(13) 


5.520(33) 


4.796(11) 










30 


13.382(33) 


10.540(33) 


8.993(29) 


7.677(39) 


7.158(35) 


6.843(27) 


6.862(23) 


31 


13.822(51) 


12.685(44) 


12.382(43) 


12.222(52) 


12.265(76) 


11.915(98) 


12.346(134) 


32 


5.743(13) 


4.781(75) 


4.070(33) 










33 


4.817(33) 


4.227(75) 


3.712(234) 










34 


8.906(13) 


8.319(39) 


8.106(19) 


8.026(29) 


8.138(27) 


8.076(23) 


8.035(19) 


35 


12.714(15) 


10.840(21) 


9.790(20) 


8.856(30) 


8.494(28) 


8.294(23) 


8.432(53) 


36 


13.604(51) 


11.073(32) 


9.597(27) 


8.365(45) 


7.762(44) 


7.499(30) 


7.293(24) 


37 


> 14.022 


13.333(100) 


12.930(64) 


12.489(98) 


12.227(186) 






38 


14.867(42) 


12.671(42) 


11.162(27) 


9.904(44)2 
9 91 3('26^3 


9 .I'W'Wl"* 


9.180(41)2 
9 153f28"l3 


8 988('30"|3 


39 


11.569(17) 


10.340(25) 


9.648(19) 


9.256(76) 


9.044(63) 


8.896(47) 


8.867(37) 


40 


>15.565 


14.487(71) 


13.733(73) 


13.034(78) 


13.001(124) 






41 


>14.189 


12.630(74) 


11.327(39) 


10.090(144) 


9.688(62) 


9.395(57) 


9.139(48) 


42 


>13.148 


12.703(64) 


12.400(41) 


11.852(96) 


11.583(106) 


11.648(178) 




43 


12.998(38) 


9.851(22) 


8.151(19) 


7.100(87) 


6.392(47) 


6.108(35) 


6.050(28) 


44 


8.363(17) 


6.836(35) 


5.920(15) 


A A6A(71 ) 


4.462(60) 


3.935(32) 


4.001(35) 


45 


7.997(13) 


6.526(21) 


5.608(9) 


4-796(81 ) 


4.211(63) 


3. 811(29) 




46 


7.960(5) 


6.534(13) 


5.684(19) 






4.608(27) 


4.351(21) 


47 


11.946(26) 


10.216(20) 


9.162(19) 


8.277(37) 


7.630(42) 


7.430(28) 


7.039(30) 


48 


11.845(28) 


10.260(29) 


9.267(24) 


8.425(38) 


7.876(49) 


7.667(32) 


7.196(28) 


49 


>13.635 


11.101(35) 


8.684(23) 


6.873(66) 




5.750(30) 


5.672(31) 


50 


13.508(24) 


12.834(27) 


12.325(26) 


11.693(64) 


11.245(102) 


11.063(85) 


11.035(164) 


51 


15.447(61) 


11.505(21) 


9.373(23) 


7.389(45) 


6.717(47) 


6.316(35) 


6.100(28) 


52 


1 94'^f25"l 






8 85n('52'l 

. C.J (J VJ ( !j 1 


8 455('37~1 


8 297(41 ~) 


7 981 ('231 


53 


>16.395 


13.645(57) 


12.169(35) 


10.855(59) 


10.433(56) 


10.163(65) 


9.886(39) 


54 


11.875(37) 


10.141(46) 


9.130(38) 


8.146(53) 


7.721(54) 


7.481(35) 


7.169(26) 


55 


14.428(29) 


12.538(25) 


11.638(24) 


10.923(57) 


10.652(69) 


10.597(77) 


10.700(72) 


56 


15.225(49) 


11.661(22) 


9.693(19) 


8.080(36) 


7.644(47) 


7.248(33) 


7.184(32) 


57 


> 14. 783 


12.960(38) 


11.574(30) 


10.204(76) 


9.551(54) 


9.270(45) 


8.958(35) 


58 


10.207(17) 


9.306(14) 


8.806(14) 


8.478(39) 


8.379(47) 


8.234(40) 


8.199(29) 


59a 


16.489(183) 


14.667(140) 


13.461(71) 


11.825(54) 


11.204(62) 


10.624(60) 


9.180(33) 


59b 


9.860(21) 


8.387(31) 


7.912(5) 


7.535(60) 


7.800(44) 


7.547(32) 


7.511(30) 


60 


8.857(13) 


8.016(19) 


7.418(15) 


6.734(30) 


6.421(38) 


6.176(22) 


5.837(16) 


61 


9.916(15) 


9.375(27) 


9.096(23) 


8.863(33) 


8.889(48) 


8.728(31) 


8.769(32) 


62a 


15.681(111) 


> 14.462 


>14.007 










62b 


15.557(66) 


>14.253 


>13.680 











^ All magnitudes from the GLIMPSE I/II/3D survey. Data in italics are from the less reliable archive deliveries. 

2 GLIMPSE I. 

3 GLIMPSE II. 
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TABLE 4 

WR Distances and Shell Sizes 



Source 


Type 


Mks 


Aks 


Av 


d 


size 






(mag) 


(mag) 


(mag) 


(kpc) 


(pc) 


9 


WN7 


-4.77 


1.27 


11.1 


9.8 


3.04 


10 


WN9h 


-5.92 


2.32 


20.3 


7.8 


1.60 


llb^ 


WC9 


-6.10^ 


1.70 


14.9 


3.8 


1.73 


16 


WN5b 


-4.77 


2.09 


18.3 


6.4 


1.89 


50 


WN6 


-4.41 


0.65 


5.7 


16.5 


1.39 


52 


WN7(h) 


-4.77 


0.66 


5.8 


5.2 


2.72 



^ Possibly not producing the shell, see text for d etails. 

2 No Mks is given for WC9 stars in Table Al of lCrowther et al.l 1(2003 ). We adopt the 
absolute magnitude of star E from their Table 9. 




Fig. 1. — MIPS images of our newly discovered circumstellar shells together with the corresponding IRAC 8fim, 2MASS Ks 

and DSS R band images. North is up and east is to the left. A 1' scale bar is indicated on the left hand side of each row of panels. 
The assumed central source for each shell is marked. The Ks and R band images are displayed at two different scales, the two right- 
most panels have a FOV of 30" X 30" . A version of the paper including the complete (online-only) Figure 1 can be downloaded from 
[http: / /web. ipac.caltech.edu/stafI/wachter/bubbles/wachter_online.pdf, 
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Fig. 2. — H and K band spectra of a sample of shell central sources indicating a late type star. The spectra have been normalized and 
offset from each other for display purpose s. For comparison the spectra of HD 207991 (K4 III) and HD 10465 (M2Ib) obtained from the 
IRTF spectral library IIHavner et al.|[2009l) are also shown. 
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2.05 2.10 2.15 2.20 

Wavelength (microns) 



Fig. 3. — J, H, and K band spectra of shell central sources characterized by H I and He I (1.70 ad 2.112 fini) emission (group 1). Star 
#10 is a newly discovered WN9h star, star #52 a newly discovered WN7 star. The spectra have been normalized and offset from each 
other for display purposes. 



18 




Fig. 4. — J, H and K band spectra of shell central sources characterized by weak or absent He I (1.70, 2.112 /^m) featur es and strong 
H I, F e II and Mg II emission (group 2A). These sources strongly resemble the stars classified as LBVs and Be/B[e] stars by IMorris et al.l 
II1996I) . The spectra have been normalized and offset from each other for display purposes. 
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Fig. 5. — J, H and K band spectra of shell central sources characterized by weak or absent He I (1.70, 2.112 /^m) features and strong Fe 
11 or Mg II emission (group 2B). The spectra are similar to those of group 2A (Figure [Jjl but exhibit weaker H I features. Note that we could 
not extract J band spectra for all of these sources. The spectra have been normalized and offset from each other for display purposes. 
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Fig. 6. — Spectra of the newly discovered WC9 (#llb) star as well as the two spectra that did not fit in the classification scheme of the 
sources in groups 1 and 2. The spectra have been normalized and ofi'set from each other for display purposes. 
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Fig. 7. — The optical spectra of the newly discovered WN type WR stars #50 (WN6o) and #52 (WN7(h)). The spectra have been 
continuum normalized and offset for display purposes. 
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6500 7000 7500 8000 

Wavelength {Angstrom) 

Fig. 8. — The spectra of the sources classified as B stars (except for source #57 which shows a F/G type spectrum). 




Fig. 9. — Color-color diagram for the sources in our sample. WR stars are indicated in red, LBVs / LBV candidates in cyan, B stars in 
blue, and F-M stars in magenta. Green symbols indicate tlie central sources that are listed in Table [T] which have not yet been observed 
spectroscopically. The comparison population of sources (black points) is composed of the photometrically most reliable point sources of a 
repres entative 1° X 1° "slice" of the Galactic plane from the GLIMPSE survey. The reddening vector was adapted from llndebetouw et al.l 



